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ABSTRACT 

The MHD stability of tangential discontinuities is first con- 
sidered. We treat these discontinuities as structured forms rather 
than as sharp breaks in the magnetic field. An unfamiliar form of 
the MHD energy principle is applied, and stability is proved pro- 
vided that there is no fluid flow tangent to the "discontinuity" plane. 
Perturbations which simply transform the system from one equi- 
librium to another are neutrally stable. Usingcomparison theorems 
we conclude that the observed stability of tangential forms in the 
solar wind implies near isotropy of the particle pressure in them. 


STABILITY OF TANGENTIAL DISCONTINUITIES 


The frequent appearance of tangential discontinuities in the solar wind flow 
(Burlaga, 1968; Siscoe et al. , 1968; Burlaga and Ness, 1969) suggests that they 
are stable structures. Contrary to the implication of the name "discontinuities," 
these structures are often thick enough that their internal details can be observed 
within the time resolution of the measurements (Siscoe et al., 1968;. Burlaga and 
Ness, 1969). They therefore are more appropriately named tangential forms . 
These forms are stationary in and convect with the solar wind. In the solar wind 
frame the magnetic field lines are straight and rotate in direction about an axis 
normal to the "discontinuity" plane. However, as stated above, the "discontinuity" 
really is thick and should be considered as a series of parallel planes (see Fig. 

1). As the field rotates through the form, it may change magnitude also. [The 
special case where the magnitude is constant is a force-free field (Woltjer, 

1958 ).] The rotation varies from monotonic, to non-uniform, to erratic (Siscoe 
et al., 1968; Burlaga and Ness, 1969). Because of this great variety of patterns, 
which have been seen to convect hundreds of earth radii (Fairfield, 1968; 

Burlaga and Ness, 1969), it seems that there should exist a universal proof of 
their stability. 

There is in fact such a proof. This proof follows from an unfamiliar form 
of the MHD energy principle. This form was published recently by Grad and 
Rebhan (1969). They show that the second order variation in total energy 
(magnetic plus plasma) of an MHD system under a trial displacement £ (£) can 


be written 


where B(r) is the unperturbed (equilibrium) magnetic field at point r, P(r) is the 
unperturbed pressure, and the integration extends over the system. The necessary 
and sufficient condition for stability is that S W > ofor all £ (r) satisfying certain 
boundary conditions. The £ (r) need not be normal modes. The energy variation 
linear in £ vanishes for equilibrium. 

'v. 

For an unbounded spatial system the energy principle exists only for the 
class of displacements £ which vanish as one goes infinitely far from the origin 
of an arbitrarily located coordinate system. Our stability analysis is therefore 
restricted to such perturbations. Practically speaking, we are thus limiting 
ourselves to perturbations of scale length short compared to the lateral (i.e., 
parallel to the discontinuity plane) extent of the observed tangential forms. 

Such is not a serious limitation, however, as tangential forms have a lateral 
extent of tens of R (Burlaga and Ness, 1969), and we expect any perturbing 
effect to be of much smaller dimension. Indeed our unbounded model is not 
even appropriate for treating perturbations of scale length comparable with the 
lateral size of the tangential forms. 

Returning to Eq. (1), we see that the first and last terms in the integrand 
are positive definite. The middle term happens to vanish for all tangential 
forms. Therefore all tangential forms are hydromagnetically stable. The 
vanishing of the middle term is easily seen from the following relations. At 
equilibrium 


2 P - syrjx b (rxt)xB --v(fj + e.?B 

V(p t f?), BVB 

Thus any situation, such as the tangential forms, where B does not vary as 
one moves along a field line, is stable. This stability is not recognizable from 
the usual (Bernstein et al., 1958) form of 8 W. 

Syrovatskii (1953) has solved the sharp boundary problem. [The solution is 
reproduced by Landau and Lifshitz (I960).] In his work, structure is ignored and 
the tangential discontinuity is represented as the sharp plane interface between 
two incompressible fluids. The magnetic field is tangent to and changes dis- 
continuous^ (in magnitude and direction) across the interface. In the absence 
of tangential fluid flows at the interface, this configuration is hydromagnetically 
stable. Our result proves stability in the more general case of an arbitrarily 
structured transition regardless of the compressibility of the fluid. 

One can not show from Eq. (1) the neutral stability of our structured system 
to perturbations which either (a) rotate the direction of B by an arbitrary angle 
uniformly in any discontinuity plane; (b) increase the magnitude of ^by an arbi- 
trary amount uniformly in any discontinuity plane; or (c) both rotate the direction 
of ^ and increase its magnitude. Any configuration obtained by processes a-c 
is another equilibrium, since for each such configuration V P = £ x g. Further, 
these equilibria are dense: there is infinitesimally close to any equilibrium 
another equilibrium obtained by rotating the magnetic field of the initial equi- 
librium through an infinitesimal angle and/or by increasing its magnitude by an 
infinitesimal amount. The transition between such infinitesimally adjacent 
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equilibria is continuous in the sense that the macroscopic variables (pressure, 
magnetic field, and density) change smoothly. One can pass by processes a-c 
between two finitely separated equilibria by a series of infinitesimal displacements 
£ through intermediate states, all of which are themselves equilibria. For each 
£, one expects 8 W = 0. 

The trouble is applying Eq. (1) to show that 8 W = 0 for such £’ s /$ ithat the ! 

~ I I 

£'s are not spatially local (in fact, £ becomes larger in magnitude as one gets 
further from the origin of coordinates) and surface terms which have been 
neglected in deriving Eq. (1) become important. 

The neutral stability of such rotational perturbations is, however, evident 
from a more primitive form (Bernstein et al., 1958) of 8W 




( 4 ) 


where 




( 5 ) 


is the force on the element of fluid (of mass density p ) displaced from its 
equilibrium position by the amounts . (In Eq. (5) y is the ratio of specific heats 
and all quantities take their equilibrium values). 

In the Appendix we discuss the form of s which produces (uniformly in each 
discontinuity plane) an infinitesimal rotation of B and infinitesimal change in|g| . 
We conclude that without loss of generality £ can be represented as the 2- 
dimensional vector 


< 
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; = x, &(i) c,(t) * 1 eu)x, - [g ?(i) t-g’d)] '[*'+ xP(i)&iij]xi I £(i) <6) 

t 

We here adopt the orthogonal coordinate system x, .Xj ,z with e t 00 = g (z)/ 

I 3 (z) | , e z along the rotation axis and therefore normal to the dificontinulty 
plane, and * e 2 = e t . P(z) and B(z) in Eq. (6) are the pressure and magnetic 
field intensity of the initial equilibrium; K' = (P' + •$') is the first order 

change in total pressure (magnetic plus thermal) produced by the perturbation 
and is independent of z (primes denote perturbations); G(z) = ft-g' B 2 -P'( yP)" 1 
is the first order change in In (|b| p" 1 ) (The density p and P are related by Pp-y 
= const.); and € (z) = (g' *e 2 ) | g I 1 is to first order the sine of the rotation 
angle. In the Appendix we also show that for £ as given by Eq. (6) F {£} and 
hence 5 W vanish. 

The importance of establishing MHD stability lies in comparison theorems 
(Kulsrud, 1964), which show that the hydromagnetic 8 W is less than any of the 
SW's obtained using one or more adiabatic invariantj in a collisionless plasma. 
The following conditions must, however, be met for the comparison theorems to 
hold: (1) the y of the gas in the MHD 8 W must bo taken as 5/3, as it has been 
in deriving Eq. (1); (2) the pressure in the adiabatic 8 W's must be taken as 
isotropic. From the comparison theorems wo thus conclude that collisionless 
tangential forms should be stable for isotropic pressures. 

If the pressure anywhere in the tangential form were sufficiently aniso- 
tropic, the mirror and firehose instabilities (Rosenbluth, 1956 j Chandrasekhar 
et al., 1958) would exist. These instabilities can be derived from the adiabatic 
SW’s and are not prohibited by the comparison theorems which hold only for 


5 



isotropic pressure. The long-term stability of the tangential forms is experi- 
mental evidence that the pressure in them is not excessively anisotropic. 
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Figure Caption 


The Structure of a Tangential Discontinuity. Five slices through the discontinuity 
are shown. B changes arbitrarily in magnitude and direction from sheet to 
sheet but always remains in the plane of a sheet. 


APPENDIX 


We here outline the derivation of the 2 -dimensional plasma displacement 
c , given by Eq. (6), which produces an infinitesimal rotation of B and an infini- 
tesimal change in ' B I . We then show that for this£ ,F {£} as given by Eq. (5) 
vanishes. We conclude by discussing an equilibrium-transforming 3-dimensional 
£ and prove that the changes in pressure and magnetic field produced by this 3- 
dimension--! plasma displacement can always be effected by a 2-dimensional £ . 

Since the perturbed state is also a tangential discontinuity we demand that 


it satisfy three conditions: 


t'j-o (A -i, 

(p+ p'jl KtK 1 (A-2) 

P t- P = (F-h p) U)- B f b ' - (E l f %) (A-3) 


(A-2) 


Condition A-l requires that the total magnetic iieid not vary as one moves along 
a field line. Condition A-2 asserts that the total pressure is everywhere constant, 
w T hiie from condition A-3 the partial pressures and § itself may change but 
only in the direction perpendicular to the discontinuity planes. In linearized 
form these conditions are 


b-vb'+ b'-pb = 0 


(A- la) 


P'+ B -S 1 - k' 


(A-2a) 




(A-3a) 






In terms of the plasma displacement £ , the pressure perturbation P ' and 

'N* 

the magnetic field perturbation B’ are (Bernstein et al., 1958) 


pL - f.v p - xP v . j 

B-7? s . VB -b?, % 


(A-4) 
(A— 5) 


We next assume that £ has only x x and x 2 components, use A-4 and A-5 
for P' and B ' ( £ «V P = £ • V g = 0 for this 2-dimensional £ ), and obtain 
expressions for the £ components from the conditions A- la - A-3a. We find 
that the most general 2-dimensional £ is 

£ = [ X-.C(i) - L Ufa,*)] $( I ) 

+ [&)X , -(&+&/[*' (A . 6) 


In (A-6) L represents a translation of plasma in the x 2 direction, uniform in 
each discontinuity plane, but varying arbitrarily from one discontinuity plane 
to another. Since a tangential discontinuity is at each z infinite and uniform in 
the x 2 direction, the translation L does not effect either a pressure or magnetic 
field change. Hence we take L = 0. Similar arguments also justify dropping H. 

When we now reintroduce A-6 into Eqs. A-4 and A-5, we find 




p= 




(A-7) 

(A-8) 
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Solving for e and G (using A- 2a for K'), we obtain 


and 


_ s'- c 


/• A* 


2 



(A- 9) 


r* 

b 


B-e' 

V 



(A-10) 


c is to first order the sine of the rotation angle. From the relation P p~ y = const., 
it is easily shown that G is the first order change in^n( | E |/c _1 ). 

Proving that F {£} as given by Eq. ( 5 ) vanishes for f as given by A-6 is a 


straightforward exercise in differentiation. The individual terms appearing in 
F are 




-O 


+ c M, P i^ PiB y fain)] 


The sum of these four terms is zero. To show this we use the facts that 
3 /3 z ( ? 1 . e 2 ) = 0 and 3 /3 z (P + B 2 /2) = 0. 

The entire procedure outlined so far in this Appendix can also be applied to 
the case where £ is assumed to have a z-component in addition to its Xj and 
x 2 components. Neglecting trivial arbitrarinesses (terms equivalent to H and 
L in A- 6), we find 
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(A-ll) 


f -X , 6 ■(i)e,(?) + fc(i)Xi- QrftB*) '[k'+xPGa) ]xz 


dc 

di 



kh) + Cd) 4 


For this f it is also straightforwardly shown that F {£} = 0. 


We note from (A-ll) that if C is independent of z so that the z-translation 
of plasma is uniform among discontinuity planes, no other component of £ is 
modified. In this case C is a trivial rigid translation of the entire system. How- 
ever, when C depends on z the plasma density and hence the plasma pressure 
are changed in a non-trivial manner by the plasma displacement C. The magnetic 
pressure is also changed since the moving, conducting plasma carries field 
lines along with it. The dependence of plasma pressure and magnetic pressure 
on density is not the same however, so that in order to maintain constant total 
pressure in the z-direction, the plasma expands (or contracts, depending on 
the sign of BC/3z) in the x 2 direction. 

The pressure and magnetic field perturbations produced by the 3-dimensional 
£ are 

'Ny 




A 

+ e, 


k'+xPQ 
HP+ B 

p'-- ~sp[s - 1-c 

L XPt 2"J 


ZlBl 


^ !g\ dtj) 

(A- 12) 


bit 

(A- 13) 


It is possible to relate e , G, and C to changes in pressure and magnetic field, 
just as was done in Eqs. A-9 and A-10 for the two dimensional £ . 
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Suppose now we define £ and G as 


€60 




lef .tf -K = /gpli^ . +. c 

~ gP-hg^ VP+B*- < 


PG 


It may then be readily verified that A-12 and A-13 become 


p / I - I Tr ~ J r XP& 

t - h IL% 2 *- e ' _ 

p'= -gf[?~ * '± i ± 2 
L XPrb- J 


Comparison with A- 7 and A- 8 indicates that these same changes B ' and P' 
can be effected by the 2 dimensional £ 

f -- x,Gc, +- ft x, - (gPt B') (k'+ XP.^JxJ Cj (A _u 

Ay *-* 


Eqs. A-14 - A- 18 illustrate the fact that the same changes in magnetic 
field and pressure can result from several different plasma motions. In 
particular changes due to an infinitesimal three dimensional plasma displace- 
ment which transforms one tangential discontinuity into another can always 
be duplicated by a ; having only Xj and x 2 components. 
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